Charge Relaxation in a Single Electron Si/SiGe Double Quantum Dot 
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We measure the interdot charge relaxation time Ti of a single electron trapped in an accumulation 
mode Si/SiGe double quantum dot. The energy level structure of the charge qubit is determined 
using photon assisted tunneling, which reveals the presence of a low lying excited state. We sys- 
tematically measure Ti as a function of detuning and interdot tunnel coupling and show that it is 
tunable over four orders of magnitude, with a maximum of 45 /is for our device configuration. 

PACS numbers: 85.35. Gv, 73.21.La, 73.63.Kv 



Semiconductor quantum dots have been widely used 
as probes of fundamental quantum physics and to im- 
plement charge and spin qubits [T^ "2 . Coherent manip- 
ulation and two-qubit entanglement have been demon- 
strated in GaAs double quantum dots (DQDs) and er- 
ror correction techniques such as dynamic decoupling 
have been employed to suppress decoherence [SUB]. As 
an alternative host material, Si holds promise for ultra- 
coherent spin qubits due to weak spin-orbit coupling, a 
centrosymmetric lattice (no piezo-phonon coupling), and 
an established route to isotopic purification [THTT]. Spin 
lifetimes of 6 seconds have been measured in Si and iso- 
topically purified ^^Si crystals can support spin coherence 
times as long as 4 seconds [71 [12] . 

While Si closely approximates a "semiconductor vac- 
uum" for electron spins, its electronic band structure 
leads to potential complications that are absent in 
the conventional GaAs/AlGaAs two-dimensional elec- 
tron gas (2DEG) system [13 . First, the ~ 3 times 
larger effective mass of electrons in Si requires deple- 
tion gate patterns to be scaled down significantly in or- 
der to achieve orbital level spacings comparable to those 
obtained in GaAs. Second, the band structure of bulk 
Si consists of six degenerate valleys, which introduces 
an additional decoherence pathway [T4|. Valley degen- 
eracy is partially lifted by uniaxial strain in a Si/SiGe 
heterostructure [15]. However, the energy splitting be- 
tween the lowest two valleys is highly sensitive to de- 
vice specifics, such as step-edges in the quantum well 
[TTJ [161 [II] • Detailed measurements of the low lying en- 
ergy level structure, and the timescales that govern en- 
ergy relaxation between these levels, are therefore needed 
in Si quantum dots [T8] . 

In this Letter, we systematically measure the interdot 
relaxation time Ti of a single electron trapped in a Si 
DQD as a function of detuning e and interdot tunnel 
coupling tc. We demonstrate a four order of magnitude 
variation in Ti using a single depletion gate and obtain Ti 
= 45 /is for weak interdot tunnel couplings [19] . We also 
use photon assisted tunneling (PAT) to probe the energy 
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Figure 1. (Color online) (a) The DQD is operated by biasing 
a global top gate at voltage Vr to accumulate carriers in the 
quantum well (left). Local depletion gates define the DQD 
confinement potential (center). Charge sensing is performed 
using rf-reflectometry (right), (b) Few electron charge stabil- 
ity diagram visible in the derivative of the reflected rf ampli- 
tude dA/dVh. (A^l,A^r) indicate the number of electrons in 
the left and right dots. (Inset) P{i,o) plotted as a function of 
detuning, for different values of Vn, showing tunable interdot 
tunnel coupling at the (1,0)-(0,1) interdot charge transition. 



level structure of the single electron system, demonstrat- 
ing spectroscopy with an energy resolution of ~ 1 /ieV. In 
contrast with single electron GaAs dots, we observe low 
lying excited states ^ 55 //eV above the ground state, an 
energy scale that is consistent with previously measured 
valley splittings [TTJ [16] . 

Measurements are performed on an accumulation 
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Figure 2. (Color online) (a) (left) Current, /, measured as 
a function of Vl and Vr near the (1,0)-(0,1) charge transi- 
tion. A cut through the finite bias triangle (right) indicates 
the presence of a low lying excited state, (b) P(i,o) plotted 
as a function of detuning e for different excitation frequencies 
/. For / ^ 15 GHz, a new PAT peak emerges (grey arrow) 
corresponding to the (1, 0)g ^ (0, l)e transition. The appear- 
ance of this PAT peak is accompanied by the suppression of 
the (l,0)g ^ (0, l)g PAT peak (black arrow) at positive de- 
tuning, (c) Transition frequencies as a function of detuning 
and (d) energy level diagram extracted from data in (c). The 
data in (c) are best fit with an interdot tunnel coupling tc = 
1.9 GHz and an excited state energy A = 55 /xeV. 



mode Si/SiGe DQD. We apply a top gate voltage Vr 
= 2 V to accumulate carriers in a Si quantum well lo- 
cated ~ 40 nm below the surface of the wafer [see Fig. 
1(a)]. The resulting 2DEG has an electron density of 
4 X lO^Vcm^ and a mobility of 70,000 cm^Vs. A 100 
nm thick layer of AI2O3 separates the top gate from the 
depletion gates, which are arranged to define a DQD and 
a single dot charge sensor. 

We first demonstrate single electron occupancy using 
radio frequency (rf ) reflectometry [20 . A single quantum 
dot is coupled to a resonant circuit with resonance fre- 
quency /r = 431.8 MHz [Fig. 1(a)] and used as a high sen- 
sitivity charge detector [21] . The reflected amplitude A is 
a sensitive function of the conductance of the single dot 
sensor, ^q, which is modulated by charge transitions in 
the DQD. We map out the DQD charge stability diagram 
in Fig. 1(b) by plotting the numerical derivative dA/dVi, 
as a function of Vi, and V^i. No charging transitions are 
observed in the lower left corner of the charge stability 
diagram, indicating that the DQD has been completely 
emptied of free electrons. We identify this charge config- 
uration as (0,0), where (A/'l,A^r) indicates the number of 
electrons in the left and right dots. 

The device is operated as a single electron charge qubit 



near the (1,0)-(0,1) interdot charge transition. Charge 
dynamics are governed by the Hamiltonian H = fcTz + 
tcCTx, where ai are the Pauli matrices. We demonstrate 
tunable interdot tunnel coupling in the single electron 
regime by measuring the left dot occupation P(i,o) ^is a 
function of detuning [Fig. 1(b), inset] [Tl[3l[22].' Qubit 
occupation is described by 
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where /cb is Boltzmann's constant, Tq ~ 100 mK is the 
electron temperature, and Q = y^e^ + 4t^ is the qubit 
energy splitting [HI |23l [24]. With V^ = 225 mV, the 
interdot charge transition is thermally broadened as 2tc 
< k^T^. Increasing t^ by adjusting V-^ leads to further 
broadening of the interdot transition. For V^ = 290, 300 
and 310 mV we extract 2tc = 3.8, 5.9 and 9.0 GHz by 
fitting the data to Eq. 1. These results show that the 
interdot tunnel coupling can be sensitively tuned in the 
single electron regime in Si. 

We investigate the DQD energy level structure in Fig. 
2(a), where we plot the current / as a function of Vi, and 
Vr with a fixed source-drain bias Vsd = 700 jaeV [25 . 
In contrast with GaAs devices, the current in the finite 
bias triangles is not a smooth function of gate voltage. 
In particular, we observe a small resonance ^ 60 fieV 
away from the interdot charge transition, suggesting the 
existence of a low-lying excited state in one of the dots. 
In a few electron GaAs DQD, orbital excited states are 
typically several meV higher in energy than the ground 
state [26 . 

Higher energy resolution is obtained using PAT spec- 
troscopy, in which microwaves drive charge transitions 
when the photon energy matches the qubit splitting, hf 
= 1], where / is the photon frequency and h is Planck's 
constant. PAT transitions are directly observed as devia- 
tions from the ground state occupation in measurements 
of P(i^o) cis a function of detuning [compare Fig. 2(b) 
and the inset to Fig. 1(b)]. For / < 15 GHz, the PAT 
peaks are symmetric around e = and shift to larger 
detuning with increasing photon energy, consistent with 
a simple two level interpretation [191 [27]. However, for 
/ > 15 GHz, an additional PAT peak emerges at nega- 
tive detuning and is not accompanied by a corresponding 
PAT peak at positive detuning. Figure 2(c) shows the ex- 
tracted transition frequencies as a function of detuning. 

The data are fit using a three level Hamiltonian that 
includes the left dot ground state (l,0)g, the right dot 
ground state (0, l)g, and a right dot excited state (0, l)e, 
as sketched in the inset of Fig. 2(b) [28 . We obtain best 
fit values of tc = 1.9 GHz and A = 55 /ieV, consistent 
with the data shown in the inset of Fig. 1(b) and Fig. 
2(a). Within the ~ 1 /ieV resolution of our measure- 
ment, we do not observe anti-crossings associated with 
(0, l)e. The energy eigenstates obtained from the model 
are plotted as a function of detuning in Fig. 2(d). For 



comparison, an excited state is observed in the left quan- 
tum dot in a second device (Device 2), with A = 64 /ieV 
[28] . For both devices, the excited state energy is highly 
sensitive to V^ and Vc^ suggesting that it is not purely 
orbital in origin [17 . 

Several additional features observed in the data are ex- 
plained by the three level model. The (0, l)g ^ (0, l)e 
intradot charge transition [dotted line. Fig. 2(c)] is not 
visible since the charge detector is only sensitive to in- 
terdot charge transitions. We also note that the (0, l)e 
^ (l,0)g PAT peak is not visible at positive detuning. 
At low temperatures, the qubit population resides in the 
ground state (0, l)g, preventing microwave transitions 
from (0,l)e to (l,0)g. Finally, the (0, l)g ^ (l,0)g PAT 
peak is suppressed when e > A due to population trap- 
ping in (0, l)e. 

We measure the interdot charge relaxation time Ti by 
applying microwaves to Vl with a 50% duty cycle and 
varying the pulse period r [Fig. 3(a)]. We focus on the 
(l,0)g ^ (0, l)g transition at negative detuning, where 
the high energy state (0, l)e is not populated. Simula- 
tions of i^(i,o) as a function of time, t, for r = 1 /is are 
shown in Fig. 3(a) for three realistic values of Ti. Dur- 
ing the first half of the pulse cycle, microwaves drive the 
(1, 0)g ^ (0, l)g charge transition, resulting in an average 
P(i 0) = 0.5. The microwave excitation is then turned off, 
leading to charge relaxation during the second half of the 
pulse cycle, with P(i,o) approaching 1 on a timescale set 
byTi. 

In the inset of Fig. 3(b), we plot P{i^o) as a function 
of detuning for r = 10 ns and r = 100 /is. As expected, 
the PAT peak is smaller for longer periods due to charge 
relaxation. Specifically, in the limit r <C Ti, there is not 
sufficient time for relaxation to occur during the second 
half of the pulse cycle, leading to a time averaged value of 
^(1,0) = 0-5. In contrast, in the limit r ^ Ti, relaxation 
happens quickly, leaving P(i,o) = 1 foi" the majority of 
the second half of the pulse cycle. Due to experimental 
limitations, such as frequency dependent attenuation in 
the coax lines and finite pulse rise times at small r, we 
are unable to drive the transitions to saturation for some 
device configurations. To extract Ti we therefore fit the 
raw i^(i,o) data as a function of r to the form 

2Ti(l _e-^/(2^i)) 

-f^(l,0) — -fmax + (^min ~ ^max) 7 (2) 

where Pmax and Pmin account for the limited visibility of 
the PAT peaks [I9l |28] . Extracted Ti values are insensi- 
tive to the rescaling of the data via Pmax and Pmin- 

The interdot charge relaxation rate is strongly depen- 
dent on the interdot tunnel coupling. This variation is 
directly visible in the data shown in Fig. 3(b) for V^ = 
225, 250 and 265 mV. To facilitate a direct comparison 
of the data, we plot the normalized electron occupation 

^(*l,0) = 0-5 + 0-25 X (P(l,0)-Pmin)/(Pmax-Pmin), USiug 

the values of Pmin and Pmax extracted from fits to Eq. 
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Figure 3. (Color online) (a) Pulse sequence used to measure 
Ti and simulated qubit response. P(i,o) = 0.5 when resonant 
microwaves drive transitions between (0, l)g and (l,0)g, and 
approaches 1 on a timescale set by Ti when the microwaves 
are turned off. (b) P(*i,o) ^^ ^ function of r extracted for 
different Vn at / = 19.5 GHz. Fits to the data give Ti = 
1.4, 0.3 and 0.05 /us for V^ = 225, 250 and 265 mV. (Inset) 
Comparison of typical PAT peaks at different r, with fixed 
Vn = 225 mV and / = 25.9 GHz. (c) Ti as a function of Vn 
in Device 1. (d) Ti as a function of Vc in Device 2. 



2 [28 . In Fig. 3(c), we plot Ti over a wide range of V^ 
for two different excitation frequencies. We see a longer 
characteristic relaxation time for larger interdot barrier 
heights, with a maximum observed value of 45 /is. The 
same overall trend is observed in data from Device 2 [Fig. 
3(d)] where the interdot tunnel coupling was tuned us- 
ing Vc' Interdot tunnel coupling is only measurable in 
charge sensing when 2tc > k^T^ [23]. For Device 1 [see 
Fig. 3(c)] we obtain 2tc = 2.4, 3.8 and 5.9 GHz for V^ = 
280, 290 and 300 mV and for Device 2 [see Fig. 3(d)] we 
obtain 2tc = 3.2 GHz for Vc = -325 mV. 

The detuning dependence of Ti is investigated in Fig. 
4(a), where we plot Ti as a function of f oc Q for the 
(0, l)g ^ {^^^)g transition [3l [22]. Data are taken at 
/ = 12.3, 19.5, 25.9 and 30.0 GHz, as indicated by the 
arrows in the energy level diagram in the upper panel 
of Fig. 4(a). Our data indicate that Ti increases weakly 
as a function of detuning for the range of frequencies 
accessible in our cryostat. 

To further investigate the excited state, we measure Ti 
for the (0, l)g -^ (l,0)g and the (0, l)e -^ (l,0)g relax- 
ation processes at the same values of / [bottom panel of 
Fig. 4(b)]. In contrast with the (0, l)g -^ (l,0)g relax- 
ation process, (0, l)e can relax via two distinct pathways 
[top panel of Fig. 4(b)]. The first relaxation process is 
a direct transition from (0, l)e -^ {^^^)g with a rate Fe, 
while the second pathway proceeds via intradot charge 
relaxation to (0, l)g with a rate Fi followed by an inter- 
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Figure 4. (Color online) (a) Ti increases weakly with / for 
the (0, l)g -^ (1, 0)g transition, (b) Ti for the (0, l)g -^ (1, 0)g 
and (0, l)e -^ (l5 0)g transitions as a function of / with Vn 
= 250 mV (lower). There are two (0, l)e -^ (l?0)g relaxation 
pathways (upper). 



dot transition to (1, 0)g with rate F*. We find that (0, l)e 
^ (1, 0)g relaxation is faster than (0, l)g -^ (1, 0)g relax- 
ation for the same energy splitting. The shorter excited 
state lifetime is consistent with either a fast direct relax- 
ation rate Fe or fast intradot relaxation followed by an 
interdot transition. Assuming Fg = Fg (since the level 
detuning is the same) and taking the measured excited 
state Ti = 55 ns at / = 21.0 GHz, we can make a rough 
lower bound estimate Fi ~ 1.5 x 10^ s~^ [29 . 

We modify the results of Raith et al. to allow the 
calculation of phonon mediated charge relaxation rates 
considering only intravalley relaxation in the far detuned 
limit {\e\ ^ tc), assuming Gaussian wavefunctions for 
the non- hybridized charge states, with dot radius a and 
dot separation 2d [28, 30j. The electron-phonon coupling 
Hamiltonian in a Si quantum well takes the form 
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Here ag^A (<^q a) ^^^ ^^^ annihilation (creation) oper- 
ators for phonons belonging to branch A (A = TAl, 
TA2 for transverse phonons and A = LA for longitudinal 
phonons) with wave vector Q, and speed of sound in Si 
cx. V is the volume of the Si quantum well layer and p 
is the density of Si. E^ and Sd are the shear and dila- 
tion deformation potential constants and Gq and Q are 
the phonon unit polarization vector and the phonon unit 
wave vector [30 . Using realistic parameters, Ti values 



calculated in this model are in order of magnitude agree- 
ment with our data [28;. However, the predictions are 
exponentially sensitive to a and d, quantities that are 
difficult to accurately determine. Moreover, the model 
predicts a relaxation rate that increases with energy split- 
ting for the range of detunings accessed in our experi- 
ment, following the power law Fi = 1/Ti cx 1]^, whereas 
we observe a rate that decreases weakly with increasing 
detuning [28]. This discrepancy may be due to a detun- 
ing dependent tc or contributions from other relaxation 
channels, such as charge noise [31 . 

In summary, we have measured charge relaxation times 
in a single electron Si/SiGe DQD, demonstrating a four 
order of magnitude variation of Ti with gate voltage. En- 
ergy level spectroscopy indicates the presence of a low- 
lying excited state. From the estimated dot radius a ~ 
38 nm, we expect orbital level spacings on the order of 1 
meV, a factor of 18 larger than the value obtained from 
PAT spectroscopy (A = 55 jaeV) [32 . This suggests that 
the low-lying excited state is a valley-orbit mixed state 
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